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ABSTRACT 

Stresses  and  velocities  arc  analyzed  for  a hollow  cylinder  of 
fuel  encased  In  metallic  cladding  and  subjected  to  high  temperature 
and  high  neutron  flux  fields.  The  material  Is  represented  by  a com- 
pressible nonlinear  thermoirradiation  viscoelastic  model.  A stress 
function  for  axisymmetr ic  plane  strain  is  Introduced,  and  the  problem 
essentially  reduces  to  solving  a nonlinear  ordinary  differential  equa- 
tion In  the  Interfacial  contact  pressure.  Some  analytical  results  are 
obtained  for  the  case  of  material  properties  Independent  of  position. 
For  the  case  ot  temperature,  flux  and  thus  material  properties 
dependent  on  position,  an  approximate  formulation  is  employed  whereby 
the  cylinder  is  divided  Into  discrete  rings  with,  constant  properties. 
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1.  INTRODUCTION 

A hollow  cylinder  of  ceramic  fuel  (e.g.,  V^O)  encased  in  metal- 
lic cladding  (e.g.,  stainless  steel)  is  a widely  used  configuration  for 
a nuclear  fueld  rod  [1].  The  cladding  provides  structural  strength 
and  also  prevents  the  fuel  from  contaminating  the  coolent,  and  there 
is  usually  a gap  between  the  fuel  and  the  cladding  at  assembly.  How- 
ever, after  the  power  of  the  reactor  has  been  on  for  some  time,  this 
gap  closes  as  a result  of  thermal  expansion  and  swelling.  An  inter- 
facial contact  pressure  is  therefore  developed  between  these  two 
materials,  and  we  may  assume  that  displacements  are  continuous  at  the 
interface.  The  plenum  pressure  is  transmitted  to  the  central  hole  of 
the  fuel,  and  this  pressure  may  be  time-dependent  due  to  the  immigra- 
tion of  gaseous  products  released  from  the  nuclear  fission  reaction. 

The  diameter  of  the  central  hole  is  assumed  to  be  constant  in  the 
present  study,  since  in  practice  its  change  is  generally  small.  The 
geometry  of  the  fuel  rod  is  shown  in  Figure  1,  where  a is  the  radius 
of  the  central  hole,  b is  the  fuel  radius,  h is  the  cladding 
thickness,  and  p(t)  is  the  plenum  pressure. 

Creep  material  parameters  are  usually  highly  dependent  on  tem- 
perature. For  problems  with  variable  temperature  fields  Morland  and 
Lee  [2]  employed  the  concept  of  a thermorheologically  simple  material, 
whereby  a temperature  dependent  thermal  reduced  time  ^ is  introduced 
to  convert  the  memory  integral  to  one  with  constant  material  properties. 
Realizing  the  analogy  between  irradiation  induced  creep  and  thermally 
induced  creep,  Cozzarelli  and  Huang  [3]  introduced  the  concept  of  a 
f luxorheologically  simple  material,  whereby  a flux  dependent  irradiation 


reduced  time  is  employed  in  the  irradiation  induced  creep  compli- 

ance function.  Also,  in  dealing  with  the  combined  case  of  thermo- 
irradiation  induced  creep  where  both  the  temperature  and  flux  fields 
are  high,  they  used  two  reduced  times  where  each  is  dependent  on  both 
temperature  and  flux.  This  approach  results  in  two  sets  of  memory 
integrals  and  differs  from  the  approach  employed  by  Rashid  [4],  where 
a single  reduced  time  dependent  on  temperature  and  flux  was  used. 

Experimental  evidence  [5]  has  shown  that  irradiation  induced  creep 
strain  for  both  ceramic  and  metallic  materials  is  essentially  linearly 
dependent  on  stress.  However,  for  thermal  creep  at  moderate  stress 
levels,  although  the  strain  for  a ceramic  material  may  still  be  linearly 
related  to  the  stress,  it  is  generally  highly  nonlinear  for  a metal. 
Therefore,  we  consider  that  the  present  problem  consists  of  a linear 
viscoelastic  cylinder  encased  in  a nonlinear  viscoelastic  cylinder. 

Because  the  fuel  rod  is  long,  this  can  be  treated  as  an  axisymmetric  plane 
strain  problem.  Also,  since  the  cladding  thickness  is  small  in  compari- 
son with  the  radius  of  the  fuel  (h  <<  b) , membrane  theory  can  be  applied 
to  the  outer  cylinder.  For  simplicity,  we  neglect  the  thermal  expansion 
and  swelling  terms  in  the  constitutive  relation,  and  also  ignore  inhomo- 
geneity and  cracking  in  the  fuel  material.  Finally,  we  note  that  we  do 
not  assume  incompressibility  since  irradiation  induced  creep  may  not 
occur  with  conservation  of  volume,  as  pointed  out  by  Gilbert  and  Straalsund 
[6]  and  Gittus  [7].  We  follow  the  approach  of  Courtine,  Cozzarelli  and 
Shaw  [8]  and  Cozzarelli  and  Huang  [3],  whereby  Poisson  coefficients  are 
introduced  for  each  component  of  creep.  For  simplicity  we  set  all  of 
these  coefficients  equal,  but  not  necessarily  equal  to  1/2  . 
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Stresses  and  velocities  In  viscoelastic  cylinders  have  been 
extensively  studied  in  connection  with  solid  rocket  tucl  propellants, 
e.g.,  see  I'l-lJ],  l'hese  analytical  studies  generally  dealt  with 
linear  incompressible  viscoelastic  materials  containing  a pressurised 
central  hole  and  in  some  cases  reinforced  by  a thin  outer  elastic  ring. 
Although  these  studies  serve  as  a useful  guide  to  the  analysis  of  the 
present  problem,  the  analytical  methods  ot  solution  employed  are  not 
directly  applicable  here  because  ot  the  Increased  complexity  ot  the 
present  fuel  rod  problem  due  t o thermal  gradients,  t lux  gradients, 
nonlinearity  and  compress  lb 1 1 1 1 y . Several  comprehensive  computet 
codes  have  been  developed  recently  tor  the  stress  analysis  ot  tuol 
rods,  employing  approximate  tormulations  and  numerical  solution  tech- 
nique such  as  finite  difference  1 1 ••  i and  finite  element  methods 
115-17],  Although  such  methods  are  necessary  tot  the  general  problem, 
the  development  ot  analytical  results  tor  various  special  cases  is 
also  Important  in  that  such  results  provide  additional  insight  into 
the  significance  ot  the  various  input  parameters.  1'he  development  ot 
"exact"  tormulations  and  analytical  solutions  is  a major  goal  ot  the 
present  paper,  but  we  do  resort  to  approximate  tormulations  and 
numerical  solutions  when  necessary. 

In  section  1,  the  basic  equilibrium,  compatibility  and  consti- 
tutive relations,  along  with  the  associated  boundary  and  initial  con- 
ditions, are  given  for  the  fuel  tod  problem.  A general  solution 
procedure  is  discussed,  whereby  a stress  tunct  ion  is  introduced  and  a 
contact  condition  at  the  fuel-cladding  interlace  is  used  to  obtain  a 


* 


nonlinear  ordinary  differential  equation  in  the  interfaclal  contact 


pressure.  Once  this  equation  is  solved,  the  stress  and  velocity  t ields 
follow  directly.  In  section  J,  this  general  procedure  is  used  to  obtain 
nond imens ional  solutions  in  three  special  cases:  1 - material  propet  ties 
independent  of  position  and  transient  creep  excluded,  11  - material 
properties  Independent  of  position  hut  with  transient  creep  included, 
and  111  - material  properties  dependent  on  position  and  transient  creep 
excluded.  In  case  1 analytical  solutions  are  obtained,  in  case  II 
"exact"  equations  tor  equal  retardation  times  are  solved  by  numerical 
techniques,  and  in  case  111  an  approximate  formulation  is  obtained  tor 
the  Incompressible  case  by  dividing  the  cylinder  into  discrete  rings 
with  constant  properties.  The  nond  imens ional  stress  and  velocity  dis- 
tributions in  the  radial  direction  and  their  variation  with  time  are 
then  discussed  for  the  three  cases  in  sect  ion  s,  with  attention  given 
to  the  effect  of  varying  the  various  nond imens ional  material  parameter 


ratios  and  the  Poisson  coefficient. 
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2 . _ BASIC  KQUAT10NS  OK  THE  PKOBLKM 


2.1  Field  Fquatioas 

The  tuel  rod  is  assumed  to  be  infinitely  long  and  in  a state  of 
axially  symmetric  plane  strain.  Thus  the  equilibrium  equation  in  polar 
coordinates  is  given  by 


r 


+ 


U) 


where  and  are  the  radial  and  transverse  stress  components 

and  r is  the  radial  coordinate.  Also,  the  compatibility  equation  is 
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where  t and  c are  the  radial  and  transverse  strain  components. 

r t? 

Finally,  the  strain-displacement  relations  are  given  as 
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where  u is  tlie  radial  displacement. 

As  discussed  in  [3],  a structural  element  under  stress  and  sub- 
jected to  high  neutron  flux  and  temperature  fields  will  experience  both 
thermal  and  irradiation  induced  creep.  A cylindrical  fuel  rod  in  a 
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reactor  is  an  element  where  both  these  effects  can  be  very  significant. 
The  total  strain  will  be  expressed  as  a summation  of  linear  elastic, 
steady  thermal  creep,  steady  irradiation  induced  creep,  transient 
thermal  creep  and  transient  irradiation  induced  creep  components.  A 
mechanical  spr ing-dashpot  model  of  the  nonlinear  generalized  Kelvin 
type  can  be  used  to  represent  such  a summation  of  strain  components 
and  this  is  shown  in  Fig.  2.  The  form  of  the  corresponding  constitutive 
relation  for  either  the  fuel  or  the  cladding  is  taken  from  [3],  and 
is  written  in  integral  form  as 
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In  order  to  limit  the  complexity  of  the  problem  we  have  assumed 
in  eq.  (5)  that  the  thermal  expansion  and  irradiation  swelling  are 
negligibly  small,  and  have  set  all  Poisson  coefficients  in  both  the 
fuel  and  cladding  equal  to  v . In  the  above  equation,  1^  = *s 


tlie  first  invariant  of  stress  and  v s.,  s.,  is  the  second 

2 2 ik  ik 
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invariant  of  the  stress  deviator  s.,  *=  a.,  - to..  6,.  . Subscripts 

lj  ij  J kk  ij  1 

Ts , Tt,  Rs  and  Rt  designate  the  thermal  steady,  thermal  transient, 
Irradiation  steady  and  irradiation  transient  creep  quantities,  respec- 


tively. The  creep  powers  and  are  non-negative 

integers,  and  the  creep  parameters  cKs » Cgt  • antl  S't  are  R08itive 

constants.  The  various  creep  compliance  functions  are  defined  as 
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where  nu(t),  n_.(t)  are  the  Irradiation  and  thermal  reduced  time 
K 1 

scales  and  tu  , t are  the  retardation  time  constants.  Also,  the 
K 1 

reduced  times  are  defined  as 
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where  4>  (T)  , Y ( <f>)  , $ (T)  and  are  some  functions  ot  tempera- 

K K 1 I 

ture  T and  flux  <f  . In  the  present  analysis  we  assume  that  T and 
have  readied  steady  state,  thus  T ” T(r)  , and  eqs.  (7. a) 


and  (7.b)  simplify  to 
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Although  the  irradiation  induced  creep  strain  will  be  taken  as 
linearly  dependent  on  the  stress,  the  thermal  creep  strain  will  In  gen- 
eral be  nonlinear.  The  constitutive  relations,  eq.  (!>)  , can  then  bo 
rewritten  lor  the  case  of  axisymmetrle  plane  strain  in  the  torm 
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where  1)  is  the  differential  operator  .l/.U  , and  where  o - v(o  +o  . ) 

z r t> 

has  been  eliminated  from  the  right  hand  side,  vurt hermore,  experimental 
evidence  [ 5 1 indicates  that,  although  thermal  creep  strain  is  nonllncat 
in  stress  lor  metals,  it  mav  be  linear  for  ceramic  materials.  Accordingly, 


in  subsequent  work  we  will  set 
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“ m for  cladding  material 
- ni^  - 0 for  core  material. 
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Following  [3],  the  axial  and  lateral  strain  for  a one-dimensional 

creep  test  o ” u U(t)  (U(t)  Is  unit  step  function)  muv  now  be  written 
o 
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where 


n " 2m+l,  X 
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It  will  be  convenient  in  later  work  to  use  n,  X,„ , \ , p_  and  p_ 

1 K i K 

instead  of  m,  CTfj,  CRs«  C . and  CRt  . 

We  consider  the  case  ol  a known  plenum  pressure  pit)  in  the 
central  hole  suddenly  applied  to  t “ 0.  The  boundary  conditions  in 
this  case  (.see  Fig.  1)  are  then  written  as 


. — . D 


10 


or(a,t)--p(t)  (1  l . *»> 

o (b,t)  - -f (t)  ( l i.bl 


where  I v. t. ) is  the  unknown  contact  pressure  at  the  intertace  ot  the 
tuel  and  the  cladding,  and  where  p(0  1 “ I (0  1 “ 0.  The  cyllndei  Is 
Initially  stress-tree  and  responds  In  a purely  elastic  manner  at  t “ 0, 
and  thus  the  Initial  condition  ntay  be  given  either  as 
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wliere  the  latter  Is  simply  the  radial  stress  of  a corresponding  elastic 
problem.  When  transient  terms  are  included  we  will  require  additional 
Initial  conditions  at  t “ 0+  ; such  conditions  may  be  generated  from 
the  governing  differential  equation  with  homogeneous  conditions  at 
t ■ 0+  and  this  will  he  discussed  In  subsequent  sections. 


2 . 2 General  Method  of  Solution 

For  the  present  plane  strain  problem  a stress  function  11  can  be 
introduced  via  the  definitions 
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thereby  ensuring  that  equilibrium  equation  fll  is  identically  sat  1st  led. 
Combining  eqs.  (2)  ,(9)  and  (.  I S > , a compatibility  equation  in  terms  ot 
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where  l’  In  the  d 1 1 t oront  In  1 operatot 
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Since  F nnd  F a to  l unci  Ions  ot  1'  nnd  which  In  tutu  nt  o 
K 1 

I unci  lottN  ot  t , dittotontlnl  equation  lib)  In  third  ordot  In  t with 

r -dependent  variable  cool  I lo I out n . The  ordot  In  time  Is  thltd  In  genet nl, 

hut  reduces  to  second  It  t - t and  to  I list  ll  t , t * (l.e., 

K l IK 

no  transient  creep).  The  solut  Ion  lot  11  with  boundary  and  lull  tal 
condlt  ions  tl  1,14)  tttav  he  ohtaitted  hv  vat  ions  analvt  leal  and  unmet  leal 
procedures,  and  this  will  he  discussed  in  detail  In  Intel  sections.  The 
stresses  then  follow  symbolically  1 row  oqs.  il5.a)  and  il‘>.h)  as 


’ - "r(t tt ) , r , t ) 


vl  .’.a) 


°,t  " 0 0 <> 1 (t)  » < • 1 ) 


t i ; .h) 


where  tit)  is  to  he  determined  hv  cons iderat (on  ol  the  cladding. 

The  cladding  is  thin  ih  " h)  and  wav  thus  he  treated  as  a mem 
hrane,  yielding  via  oq.  tll.h) 
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°r  » -f(t) 


o-  “ r f(t) 


e h 


(18. a) 

(18. b) 


where  superscript  c designates  a variable  in  the  cladding.  (A 

variable  without  a superscript  will  always  relate  to  the  fuel.)  Now, 

•or  plain  strain  the  invariant  tern,  in  eqs.  (9. a)  and  (9.b)  simplifies 
to 


j + 1-2v  ,2  1 

2 6(l+y)  1*2 
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(19) 


Substituting  <18..)  <18.b)  ,„to  eq.  „0  obtaJn  for  (he 

cladding 
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Then,  substituting  this  result  into  constitutive  relation  (9.b),  we 
obtain  the  relation  for  the  tangential  strain  in  the  cladding  as 
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with  modified  constants 
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are  assumed  Independent  ot  r since  the  cladding 


where  and  F^, 

is  modeled  as  a membrane.  Finally,  we  assume  that  the  tuel  and  the 
cladding  press  against  each  other  tightly  such  that  the  radial  dis- 
placement at  the  Interlace  is  continuous,  and  thus  we  obtain  t rom 
eq.  (4)  the  contact  condition 


The  formulation  ot  the  problem  is  now  complete,  since  a suhstltu 
t ion  of  eqs.  (9 .a) , (9 .b)  , ( 17  .a)  and  (17, b)  into  eq.  llb.af  yields  a 
nonlinear  differential  equation  in  t . Once  this  equal  ion  Is  solved, 
stresses  follow  directly  t rom  eqs.  fl/.a.b)  and  the  radial  velocity  mav 
be  obtained  from  eq.  i.4)  in  conjunction  with  eq.  This  general 

procedure  will  be  applied  in  the  next  section  to  the  lollowing  thtee 
special  cases  ot  tuel  rod  problems:  1 - material  properties  indepen- 
dent ot  position  with  no  transient  creep  components,  ll  - material 
properties  independent  ot  posit  ion  but  w i t h transient  creep  components, 
and  ill  - material  properties  dependent  on  position  with  no  tianslent 


creep  components. 


J.  SOLUTIONS  IN  SOME  SPECIAL  CASES 


As  previously  mentioned,  the  thermal  and  irradiation  creep  com- 
pliance functions  are  functions  of  temperature  T and  neutron  flux  if  , 
which  in  turn  are  generally  dependent  on  the  position  r . In  order 

to  account  for  this  dependence,  we  express  the  functions  F_  and 

R 

E^.  [see  eqs.  (8)]  in  the  form  of  M-th  order  radial  power  series 
expansions  about  their  mean  values,  i.e., 

FR(r)  ‘ FRorR(r)  = FRo  lQ  \i  ^ (23‘a> 

fT(r)  *=  1'ToIT(r)  “ l kTi  r (23.b) 

i-0 

where  FRq  and  F^o  are  the  mean  values  of  FR(r)  and  F^(r)  in 

the  fuel,  which  in  turn  require  that  the  mean  values  of  F and  T_ 

K 1 

equal  unity.  In  the  previously  described  Cases  I and  11  F » T = 1 

R T 

and  thus  r„  - F and  F * F , but  in  Case  III  F„  and  F„  are 
K KO  1 To  R I 

not  constant  but  are  taken  as  functions  of  r in  accordance  with 
eqs . (23) . 

3.1  Case  I - Material  Properties  Independent  of  Position  - Transient 
Creep  Terms  Excluded 

In  this  case  m^  - m^t  *■  0 and  tt  • tr  “ » and  constitutive 
relations  (9)  simplify  for  the  fuel  to 
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Accordinglv,  the  compatibility  equation  in  the  stress  function  H 
[eq.  (lb)]  reduces  to 
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This  equation  may  be  solved  by  successive  integration  to  yield 


H c c.(t)  7—  + c~(t)inr  + c,(t)+  c,(,r)e  u,'t  (2b) 

(CRs+  CTs)E 

where  u>  •>  and  c^(t),  c.,  (t)  , c^(t)  and  c,  are  func- 

tions of  integration. 

Using  eq.  (2o)  in  definition  (15. a)  with  a prime  indicating  d/dr 
we  obtain  for  the  radial  stress 
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^,(t)  c (t) 

— x + — + 


from  which  we  may  conclude  that  c,fO~)  - e,(0“)  - c’fr)  - 0 in 
order  that  initial  condition  (14. a)  be  satisfied  at  every  r . We 
may  now  enforce  boundary  condition  (13)  with  pfO-)  - £ (0“)  - 0 
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to  evaluate  integration  functions  c^t)  and  c2(t) . and  the  stresses 
are  obtained  with  the  use  ot  eqs.  (15)  as 
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Note  that  eqs.  (17)  are  In  the  same  torm  as  the  elastic  solutions 
[e.g.,  see  Kq.  (14.b)l,  but  here  f(t)  is  strongly  dependent  on  the 
creep  properties.  The  circumferential  strain  at  r - b in  the  fuel 
now  tollovs  from  eqs.  (14. b)  and  (17)  as 


fl+V 
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(28) 

Equating  eqs.  U8)  and  (11),  we  finally  obtain  for  the  contact  pres- 
sure the  nonlinear  ordinary  differential  equation 
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(29. a) 


where  a dot  Indicates  d/dt  and 
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1 
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K ib- -a- 1 

-a-  + C ) (1-v) 
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(b--a-)(l+v) 


t,J9.b-t  ) 


with  subscript  1 used  to  identity  quantities  in 
The  pressure  in  the  central  hole  of  the  tuel 
constant  after  being  suddenly  applied  to  t - 0 , 
In  order  to  obtain  the  solut  ion  to  eq.  (,J9.a)  for 
initial  value  of  f(t)  . This  can  be  obtained  by 
with  respect  to  t from  0-  to  0+  and  setting 
the  problem  is  reformulated  as 


Case  1 . 

is  now  assumed  as 

i.e. , pit)  - p l it) . 

o 

fit),  we  require  an 
integrating  eq.  (29) 
f(0~)  " 0 . Then , 


f + Cj  f + K f 
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with  the  initial  condition 
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It  is  difficult  to  assign  typical  numerical  values  to  the  coeffi- 
cients in  eq.  fJO) , since  they  contain  creep  constants  which  vary  ovet 
very  large  ranges.  However,  it  this  equation  is  rewritten  in  nondimen- 
sional  form  we  need  only  specify  ratios  ot  these  creep  constants,  and 
the  analysis  of  the  solution  is  greatly  facilitated.  Thus  we  Introduce 
the  nondlmensional  variables 
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and  eqs.  (30,31)  become 
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where  the  dot  now  Indicates  d dt  . In  eq.  (33)  the  nondimens lonal 
coefficients  are  given  by 
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are  nond Imensional  geometric  constants;  and 
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are  nondimens ional  material  constants.  Note  that  and  K,  have 

been  expressed  in  terms  of  one-u  nnens ional  creep  constants  through  the 
use  of  eqs.  (12.a-d)  in  conjunction  with  definitions  (21. a,  24. c). 

The  solution  to  differential  equation  (33)  is  easily  obtained  in 
implicit  form  as 
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Since  f -*  0 as  t -*•  05  , an  algebraic  equation  in  the  asymptotic  con- 
tact pressure  follows  from  eq.  (33)  as 


f («>) 


_a] 


V(° 


.)2"1 


0,1,2, 


(40) 


In  the  special  linear  case  m = 0 an  explicit  solution  is  obtained 
from  eq.  (33)  as 
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where  we  see  that 
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We  may  now  express  the  nondimensional  stress  and  velocity  solutions 


in  terms  of  f . First  we  let 
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where  r,  o , oq  and  u are  the  nondimens ional  radial  coordinate, 
radial  stress,  tangential  stress  and  radial  displacement  respectively, 
with  f(0+)  given  in  eq.  (31).  Eqs.  (27)  then  yield 
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where  we  have  introduced  the  nondimens ional  constant 


v = (1-v+vh) (l-a2)+(a2+l-2v)Eh  .... 

2(l-v)a2EH  C ; 

The  initial  values  of  stress  o (r,0+),  o„(r,0+)  and  the  asymptotic 

r v 

values  a (r,°°),  a.(r,«>)  follow  directly  from  eqs.  (34), (40)  and  (43). 
r a 

Finally,  eqs.  (4),(24.b)  and  (23)  give  for  the  nondimensional  velocity 
v = 3u/3t 


a2(l+v)[,  „ V-,  ll  n-  , ?2m+ri  f(l-2v)r  , fl) 

V = (1-~S~2)  \k  (L"2^)r+  = " l+(l-aI)f-6If  I — p + - j( 45) 


where  the  initial  and  asymptotic  values  v(r,0+)  and  v(r,“>)  again 
follow  from  eqs.  (34)  and  (40)  respectively. 

In  the  next  subsection  we  investigate  the  effect  of  the  transient 
creep  terms. 
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Altov  using  Initial  conditions  lUn,48l  to  show  that  . it'*  * . o'*  - 0 
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in  the  cladding  equal  to  those  in  fuel  (i.e.,  t = - ?£) , we 

1 K 

obtain  from  eqs.  (21,50)  the  second  order  nonlinear  differential  equa- 
tion 
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with  subscript  11  indicating  constants  in  Case  11. 

In  order  to  find  the  solution  for  f , we  require  two  initial 
conditions,  which  are  readily  obtained  by  integrating  eq.  (51)  twice 
with  respect  to  t from  0~  to  0+  and  then  setting  t'(0_)  » f(0“)-0. 
it  the  pressure  at  the  center  hole  Is  again  maintained  constant  at 
P for  t N 0 , we  obtain  the  following  differential  equation  and 
initial  conditions: 
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(52) 

151) 
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Note  that  ^jj/A  ■ H ^ / B ^ [see  eqs.  (51)  and  (29)  j,  and  thus  the 
initial  contact  pressure  given  by  eq.  (53)  is  Identical  with  the  result 
given  by  eq.  (31)  in  Case  I.  A closed-form  solution  to  eq.  (52)  Is 
not  generally  obtainable,  but  after  nondlmenslonallzat ion  a numerical 
solution  may  be  obtained. 

We  shall  employ  the  same  definitions  for  the  nond imetis iona 1 contact 
pressure  and  the  nondimenslonal  time  as  in  Case  1 [eqs.  (32.a,b)J.  Eqs. 
(52-54)  then  become 
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Note  that  Case  11  degenerates  to  Case  1 as  f ->  « , i.e.,  for 
K ^ ■ 0 eq . (55)  reduces  to  eq.  (33).  A standard  com- 

puter subroutine  was  used  at  the  SUNYAB  Computer  Center  to  obtain 
numerical  solutionsto  nonlinear  differential  equation  (55).  The  routine 
uses  a predictor-corrector  method  for  nonstiff  equations,  and  a 
variable-order  method  for  stiff  equations.  The  latter  method  is  used 
in  dealing  with  our  problem  when  m is  large.  The  asymptotic  value 
?(*»)  may  be  obtained  analytically  by  setting  F - f » 0 in  eq.  (55), 

* “ ?m+l 

whereupon  we  get  a^f^’)  + 6 f (oo) " ■»  which  can  be  shown  to  be 

identical  with  eq.  (40)  of  Case  1. 

We  now  introduce  the  same  definitions  for  the  nondimensional 

variables  r,  o , o and  u as  in  Case  1 Isee  eqs.  (42.a-d)],  with  the 
r y 

additional  definition  for  the  nondimensional  strain 


(t>4) 


Using  these  definitions  we  of  course  obtain  from  eqs.  (27)  the  same 
form  for  the  nondimensional  stress  relations  as  in  Case  1 li.e.,  eqs. 
(43)].  The  governing  differential  equation  for  the  nondimensional 
circumferential  strain  is  now  obtained  by  substituting  definition  (b4) 
and  stress  solutions  (43)  into  eq.  (4b. b),  and  thus  we  obtain 
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3 3e„  -2 


3t 


8 + K “"8  _ a2 (1+y) 


33E  (1- 


+ frji-wa^rt!! 


+ K3(1-2v+  -=7)k 


G(r,t)  (65) 


where  G(r,t)  is  known  in  terms  of  the  numerical  solution  for  f(t). 
The  initial  conditions  for  e ore  obtained  by  following  the  same 

U 

procedure  used  in  obtaining  f(0+)  and  f(0+)  for  eq.  (52),  and 
accordingly  we  obtain 


/-  r,+N  (1+v) 

i(r’°  > - O^T 


k(l-2v+  ^)~(.^=*r  + fT) 


(66) 


3 - ^+s_  (1+v) a2 

at  e6(r*°  )=  ll^y 


(en  + Yir2"K5)(1i^i  + 


+k(l+K5)  (l-2v  + 4-r) 


(67) 


Finally,  eq.  (65)  coupled  with  condition  (67)  can  be  solved  by  the 

method  of  integrating  factor  to  obtain  3e  /8t  . By  the  use  of  eq.  (4), 

u 

we  then  find  the  nondimens ional  velocity  as 


v 


3u 

3t 


-£e 

r af 


-K.t 

lt'e(;-0+)e  + 


-K  (t-t') 
e J G(r 


,t’)dt’J 


(68) 


Note  that  the  initial  value  of  v(r,0+)  obtained  with  the  use  of 
eq.  (67)  is  not  the  same  as  the  result  obtained  by  setting  t = 0 
in  eq.  (45)  of  Case  I.  However,  one  can  show  that  the  fesymptotic 
value  v(r,<»)  [obtained  from  eq.  (65)  by  setting  32£g/at2  = 0 , 

ai./at  = v(r,<*)/r,  f = f » 0 and  f = f(°°)]  is  identical  with  the 

0 

result  obtained  in  Case  I by  letting  t ■+•  00  in  eq.  (45) . 


— ! — - 
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All  solutions  obtained  in  tills  section  reduce  to  the  corresponding 

solutions  for  Case  1,  when  we  let  f ♦ ® and  set  K.  “ K,  - K,  - K “ 0 . 

i <*  j b 

In  the  following  subsection  we  again  ignore  transient  creep,  but  intro- 
duce the  eftect  of  position  dependent  material  properties. 

3.3  Case  III  - Material  Properties  Impendent  on  Position  - Transient 
Creep  Terms  Excluded 

In  this  case,  the  material  properties  are  not  constant  but  are 
taken  as  dependent  on  the  position  in  accordance  with  power  series  (23). 
Tlie  constitutive  relations  for  the  fuel  are  thus  obtained  from  eq.  (d) 
cl  s 


~ + ^][(l-v)or  - voe] 


(bd.al 


e6  ■ pr  + ZDEi]t(1’v>0e  " V°r 1 


(69. b) 


where  using  eq . (24.c) 


K(r)  - C^yr)  * eTsrT(r) 


(701 


which  is  a prescribed  function  of  r . In  the  present  case,  compati- 
bility condition  (16)  becomes 


3 f j./  Dr,,  \ 3"H  v 3H 

W[\n  + Mr)J  La'v)7T  ■ r 


VH  1 3H~|1 


3r‘  r 3r, 


(71. a> 
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where 


A;  f ) 


mil 

1+v 


(71. b) 


The  general  solution  to  eq.  (71. a)  is  not  easily  obtained,  but 
in  the  special  case  of  incompressible  materials  (v“l/2)  it  simplifies 
to  the  factored  form 


(72) 


Eq.  (72)  is  now  readily  solved  through  the  use  of  an  integrating  factor 
and  successive  integration,  and  thus  we  confine  our  attention  to  this 
special  case  in  the  remainder  of  this  subsection.  Accordingly,  the 
stresses  as  defined  in  eqs.  (15)  are  obtained  as 

r . f -A(r ' ) (t-t ' ) 

°r  " / 7^  i g2(t')e  dt'dr' 

a 0+ 

f -A(r ' ) t 

+ j g4(r')e  dr'  + gj(t)  l7J.a) 


r . r -A(r ') (t-t ')  f -A(r ' ) t 

6 ” J j g2(t')e  dt'dr'  + 84(r’)e  dr'  + g^t) 

a 0+ 

. ) -A(r ' ) ( t-t ' ) -Afr)t 

+ j g,(t’)e  dt'  + g4(r)e  (73. b) 

0+ 


where  gjU),  g,(t)  and  gA(r)  are  three  integration  functions. 
Expressions  for  these  integration  functions  may  be  obtained  by  applying 


JO 


conditions  (.13)  and  ils.b'  to  eq.  (73. a),  and  accordingly 


SjU)  - -pit) 


K4  <• r > “ §7~a? I P (0+)  - 1 (,0+)  K~r) 


b t 

f JL_  f „ 

J r'  ' bJt,t  )e 
a 0+ 


-A<r')(t-t') 


dt'dr'  - pit)  — t it)—  p7p-[p(0+)-i  10+) ) 


2 -Ur'H 

«** 


l'e  seo  lhat  t hough  g^t)  and  g^ir)  are  now  determined  via  eqs. 
J -.-75),  eq.  l/o)  must  tirst  be  combined  with  the  contact  condition 
before  g.u'*  and  tit)  can  be  determined.  The  asymptotic  values 
°r ^ 1 ’*)  Jllj  ^U'.*)  are  easily  determined  by  setting  3/Jt  - 0 
m eq.  i 2)  and  repeating  the  above  solution  procedure,  yielding 


o U'  »‘v) 
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■ sa*’) 

J r^s(v')  dr'  + *!<*> 


n>(") 

) tr,*>)  + —5 

r r-.\(r) 


f77.a> 


1 77 . b) 


In  the  present  case  the  transient  creep  terms  are  absent  and 
the  Poisson  coefficients  are  all  equal  to  1/2  ; the  equation  tor  the 
c ircumt eretuial  strain  in  the  cladding  then  follows  from  eq.  (dl)  as 


»*  - [*  -3  . Ks  fa  b.Jm.Jm  kf 

*9  1P’C  dT  + p^<1  + K>  1 /1+K>f 


Also,  the  circumferential  strain  in  the  fuel  at  r - b is  obtained 


from  eq.  tb9.b),  which  becomes  for  v » 1/2 
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\>(r) 


u 


o -o 
t»  r. 
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Then  by  substituting  eqs.  (73)  tor  tin*  stresses  into  eq.  f/9),  we 
Set  with  the  use  ot  definition  (71. b)  the  simple  result 


£e(r) 


j 

4E?  l1 


v t > 


tBO' 


wtiere  we  se  • that  the  assumption  ot  incompressibility  leads  to  a 
strain-rate  which  is  separable  in  r and  t . By  equating  eqs.  (7$) 
and  (BO)  at  r - b tin  accordance  with  contact  condition  fJd))  we  are 
able  to  obtain  the  expression 


B,U)  " 


Eb 

3 


'■'L'1*  bi + 


Is  _ b ' 
,Jm+lU+  h^ 


Jm+ 1 


tBl  ^ 


which  may  now  be  combined  with  eq.  t7c'>. 

Because  the  integrand  on  the  lett  hand  side  ot  eq.  to'  is  not 
separable  in  r and  t , repeated  different lat ion  with  respect  to  t 
will  not  yield  a simple  differential  relation  between  g.tt'  and 
ttt).  In  order  to  obtain  such  a differential  relation  which  in  turn 
would  yield  a differential  equation  for  the  contact  pressure  when  com- 
bined with  eq.  tBl)  , we  employ  an  approximate  formulation  in  which 
the  fuel  core  is  divided  into  N concentric  rings.  In  each  ting  the 
material  properties  are  assumed  to  be  constant,  and  we  set 


•Hr)  - A , r . > r r . (i  - l.'-’.N;  r - a.r  - b' 
1 i-l  1 ON 


where  A^ 


is  the  mean  value  in  the  i-th  ring.  Using  eq.  tBA'  we 


may  approximate  eq.  ( ~t>)  .is 


pit)  - lit)  - 


l S At) 

1-1 
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where 


' 1 , 1 -Aj(t-t ')  -A,t 

Stu)  - r-n  »2<f>*  dt  'dr '+~e  R, 


r 0+ 

i-1 


lp(0+)-f  (0+> 


- I dr'  11  1 , 

i I r 2 r,  r- 


‘ ri  rl-l 


1.8  1 . b-c  3 


The  integrals  in  the  functions  are  now  separated  in  r and  t . 

and  thus  the  K-th  order  time  derivatives  of  eq  • (8.3. a)  with  1 v K «■  N 

a ssume  the  si mp 1 e to rm 


V Aj.<  + i ' I k- 

i-1  J = 1 ' ~(«1  -I 


1 , • • • ,N 


where  superscript  (. K)  indicates  the  K-th  time  derivative  ot  a vat  lab  > 
Kqs.  (.83. a)  and  (8s)  constitute  a set  ot  N+l  equations  in  the  N 
functions  S^(t)  and  the  Integration  function  g.yt).  This  set  mav 
be  reduced  by  simple  algebraic  elimination  to  a single  ditteientt.il 
relation  between  g,(t)  and  t ( t > . with  all  the  S.yt)  eliminated. 

As  a simple  illustration  we  consider  a two-ring  approximation  y i . e . . 
N-2)  with  rings  of  equal  thickness;  the  three  equations  given  by  eqs. 


(83. a)  and  (84) are  now 
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P ~ f 


l Si 

i-1 


— — 

p _ f - - £ A1Si  + g2(t)  l Ri 


i-1 


i-1 


p - f 


l A*S  - g (t)  l A R + g2(t)  l Ri 
i-1  i-1  i*=l 


(85) 


(86) 


(87) 


After  eliminating  and  S0  from  the  above  equations,  we  obtain 

the  following  differential  relation  between  g,,  and  f : 

(p-f)+(A1+A2)(p-f)+A1A2(p-f)»(R1+R2)g2(t)+(A1R2+A2R1)g2(t)  (88) 

This  relation  may  now  be  combined  with  contact  condition  (81)  to  obtain 
the  governing  differential  equation  in  the  contact  pressure  for  the  two- 
ring  approximation.  Accordingly,  we  obtain 


,.2m+l.  2m+l 

+ Bin£  + cm£  + Dm(£  )+  Em£ 


Snp  + hiiip  + ImI> 


(89. a) 


where 


Am  “ (Ri+  R2)b'  b (1+  b + 1 
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<-Wb2  7^i<1+  r> 

e ; 

(AlR2+A2Rl>b2  77^T(1+ 

3 • 2 


Gm  - 1 


HIII  = (A1+V 


IIII  = A1  A2 


(89 .b-i) 


with  subscript  III  indicating  constants  in  Case  III. 

Since  eq.  (89. a)  is  second  order  we  require  two  initial  condi- 
tions, which  are  obtained  by  the  same  procedure  discussed  for  the 
previous  cases  as 


f(0+)  = 


f(0+)  = 


GIIIPo 


HIIIPo 


BIIIGIIIP0 


DIIl/GIIIP0 

AIII 


(90. a) 


(90. b) 


For  additional  simplicity  we  shall  also  assume  that  the  functions 


TR(r)  and  ^(r)  in  eq.  (70)  are  equal  in  each  ring  of  the  two-ring 
approximation.  Then  in  eqs.  (89. b-i)  the  mean  value  of  A(r)  in  the 
i-th  ring  becomes  upon  setting  v = 1/2  in  eq.  (71. b) 


Ai  " Fi 


(CRs  + £Is)2E 


i-1,2 


where  T.  is  the  corresponding  mean  value  of  r (r)  (or  T (r)) 

1 K i 


in  the  i-th  ring. 
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For  the  load  we  assume  again  that  a constant  pressure  P0  Is  maintained 
at  td»e  central  hole  for  t>0.  After  introducing  the  same  nondimens ional 
contact  pressure  f and  time  t as  in  Case  I [eqs.  (32)],  eqs . (89. a) 
and  (90)  became 
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In  the  above  ( ) = d/dt  , and 
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For  I’  « I'  " 1 second  order  differential  equation  (92. a)  is 
equivalent  to  first  order  equation  (33)  of  Case  1 (with  v - 1/2)  , 
since  it  may  be  show  that  then  eq.  (92. a)  is  a linear  combination  of 
eq.  (33)  and  its  first  derivative.  We  also  see  that  eq . (92. a)  is 
very  similar  in  iorm  to  nonlinear  second  order  equation  (55)  of  Case  11, 
and  thus  the  same  numerical  procedure  may  be  used  in  obtaining  F . 

Had  we  taken  more  than  two  rings  (i.e.,  N ' 2)  in  our  approximate 
formulation,  we  would  have  obtained  an  N-lh  order  differential  equation 
in  t . However,  the  same  numerical  method  can  also  be  applied  to  the 
solution  ot  this  higher  order  equation.  Additional  initial  conditions 
must  ot  course  be  generated  in  accordance  with  the  procedure  pre- 
viously described.  As  in  the  previous  cases  the  asymptotic  value  f(>v) 
may  be  obtained  analytically  by  setting  I « f * 0 in  eq . (92. a), 
yielding  ^ | [ j + 1 1 l * ' (“')  ” | i i which  may  then  be  used  in 

eqs.  (77).  Note  that  since  the  present  material  is  not  equivalent  to 
one  with  constant  properties  as  t » , the  asymptotic  value  will 

in  general  be  different  from  those  in  the  previous  cases. 

Having  obtained  approximate  values  for  the  contact  pressure  and 

the  integrating  function  g,(t)  via  oq.  (81),  the  stresses  are  then 

obtained  from  eqs.  (73-75).  The  best  procedure  would  be  to  drop  the 

multi-ring  approximation  at  this  point,  and  evaluate  the  integrals  in 

eqs.  (73)  numerically  over  r using  the  actual  distribution  for  A(r) 

and  thereby  obtain  continuous  distributions  for  both  o (r)  and  o,(r). 

r 0 

However,  a simpler  estimate  to  o may  be  obtained  by  also  employing 
the  step  approximation  t o A(r)  in  eq . (73. a).  Thus  we  nond imens ion- 
all.-e  in  accordance  with  eqs.  (42)  and  employ  the  two-ring  approximations 


to  g,(t)  and  A(r)  to  obtain  the  est  imate 
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Note  that  a is  continuous  at  the  junction  of  the  regions,  whereas 
the  slope  do  /'Jr  is  discontinuous.  A similar  approximate  procedure 
applied  to  eq.  (73. b)  lor  a would  yield  a less  useful  result, 
since  we  see  from  the  last  term  that  oQ  would  be  discontinuous  at 
the  junction  of  the  rings. 

Turning  now  to  the  evaluation  of  displacement  , we  first  com- 
bine eqs.  (80)  and  (81)  to  obtain  the  strain-rate  in  terms  of  f as 
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After  nondimens ionalizing  in  accordance  with  eqs.  (32)  and  (64)  this 
relation  becomes 
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Note  that  the  tangential  strain-rate  and  its  derivatives  are  continu- 
ous throughout  the  whole  fuel  region.  Finally,  the  nondimens ional  radial 
velocity  is  obtained  from  eq.  (4)  as 


r 


at 


(.97) 


where  De^/at  has  been  given  in  eq.  (9b).  The  initial  and  asymptotic 
values  v(r,0+)  and  v(r,«0 


follow  in  the  usual  manner. 
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In  the  next  section,  some  typical  values  are  chosen  for  the 
nondimensional  material  parameter  ratios  and  the  Poisson  coefficient, 
in  order  to  study  the  nature  of  the  stress  and  velocity  distributions 
in  the  three  parts  of  this  section. 


4.  DISCUSSION  OF  RESULTS 


In  order  to  illustrate  the  solutions  obtained  in  the  previously 
discussed  Cases  1,  II  and  III,  we  select  for  the  nond imensional  geo- 
metric and  elastic  constants  the  typical  values  a ■ 0.3,  h “ 0.1 
and  E * 1.0  [see  eqs.  (37. a, b)  and  (38.a)].  Turning  first  to 
Case  1 (material  properties  independent  of  position  - transient  creep 
excluded)  , we  shall  select  two  sets  of  values  for  the  nond imensional 
creep  constants  and  K.,  in  recognition  of  the  large  range  of 

variation  for  these  constants.  Using  the  one-dimensional  test  con- 
ducted at  a stress  equal  to  the  internal  pressure  pQ  as  the  measure 
of  material  behavior,  we  shall  designate  as  Case  I. a the  situation 
in  which  the  steady  thermal  creep  strain  rates  are  of  the  same  order 
in  the  fuel  and  cladding  materials,  while  being  of  one  order  higher 
than  the  corresponding  steady  irradiation  creep  strain  rates.  Thus, 
referring  to  eq.  (11. a),  we  set  in  Case  l.a 
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whereupon  eqs.  (38.b,c)  yield  the  values  *=  0.091  and  K,,  ■=  0.909. 
In  addition,  we  consider  the  case  in  which  the  steady  and  irradiation 
creep  strain  rates  are  of  the  same  order  in  the  fuel  material , and 
of  one  order  higher  than  the  corresponding  creep  rates  in  the  cladding 
material.  Designating  this  as  Case  I.b,  we  accordingly  set 


rT. . "ro  . "T.C1 
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0.05. 


and  obtain  - K,  - 0.05. 

The  nondimens ional  stress  and  velocity  solutions  in  the  fuel  for 
Case  I. a are  given  by  the  solid  curves  in  Figs.  3-to.  For  these  and 
all  subsequent  curves  0.3  <_  r 1.0  (from  central  hole  to  interface) 
and  0+  <.  t < " (from  initial  to  asymptotic  solutions),  for  m « 0 
and  2 (linear  and  nonlinear  thermal  creep  in  cladding)  and  with  v » 0 
and  0.5  (non-contracting  and  incompressible  materials).  Case  l.a  is 
especially  interesting  'n  that  for  m ■ 0 we  obtain  ■ 6^  ■ 1 

as  a consequence  of  + K,  - 1 [see  eqs.  (35-3b)],  and  thus  eq. 

(41. a)  gives  a constant  contact  pressure  f - 1 . Fig.  3 gives  the 
radial  and  circumferential  stress  distributions  plotted  versus 
position  r in  the  fuel;  the  value  of  v had  a minor  effect  on  the 

stresses  in  tills  case  and  thus  only  the  curves  for  v “ 0.5  have 

been  shown  here.  We  see  that  is  compressive  and  is  tensile 

at  all  points.  Also  we  find  that  the  initial  elastic  and  asymptotic 
stress  solutions  for  m ■ 0 are  identical,  which  is  a direct  con- 
sequence of  the  fact  that  f ■ 1 . Figures  4 and  5 present  the  radial 
and  circumferential  stresses  respectively  at  the  interfacial  position 
r ■ 1 , but  plotted  now  versus  time  t . Note  that  for  m j4  0 
la  | " f increases  with  time  whereas  cr.  decreases  with  time,  and 

i r 1 e 

that  the  effect  of  v is  minor  when  compared  to  the  effect  of  m . 
Finally,  the  radial  velocity  at  the  interface  is  plotted  versus  time 
in  Fig.  b,  and  we  see  that  the  values  of  both  v and  m have  a very 
significant  effect  on  these  time  decreasing  solutions. 

A corresponding  set  of  stress  and  velocity  solutions  for  Case  l.b 
is  given  in  Figs.  7-10.  This  case  differs  sharply  from  the  previous 
case  in  that  f always  grows  with  time  (see  Fig.  8)  and  is  never 
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constant.  Furthermore,  Figs.  7-9  show  that  the  value  of  v now  has 
a very  significant  effect  on  the  stress  solutions.  In  fact,  the 
Poisson  ratio  plays  a very  interesting  roie  in  Case  I.b  in  that  for 
v = 0.1822  the  solutions  are  independent  of  the  creep  power  m , and 
curves  for  this  value  of  v have  also  been  shown.  On  comparing 
Cases  I. a and  I.b  we  find  that  the  various  solutions  approach  differ- 
ent asymptotic  values,  and  the  rate  of  approach  is  generally  faster  in 
Case  I. a.  We  also  note  that  although  the  initial  stresses  are  identical 
in  both  cases  the  initial  velocities  are  different  since  these  are 
affected  by  the  stress  rate. 

We  turn  next  to  Case  II  in  which  the  material  properties  are 
independent  of  position  while  transient  creep  terms  with  equal  retarda- 
tion times  are  included,  and  for  the  sake  of  brevity  we  will  consider 
only  one  set  of  material  creep  constants.  We  again  use  the  one- 
dimensional test  at  p as  the  measure  of  material  behavior,  but  now 

o 

also  specify  that  the  creep  terms  in  eq.  (11. a)  be  compared  at  t = 2t 
since  the  transient  terms  have  essentially  reached  their  asymptotic 
values  at  this  time.  We  shall  take  Case  II  as  an  extension  of  Case  I. a, 
and  thus  we  assume  that  under  the  above  conditions  all  of  the  thermal 
creep  strain  components  in  the  fuel  and  cladding  materials  are  of  the 
same  order  as  the  elastic  strain,  whereas  all  of  the  irradiation  creep 
strain  components  are  of  one  order  lower.  Accordingly,  referring  to 
eq . (11. a)  we  set 
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whereupon  eqs.  (38.b,c)  aiul  (bl.a-d)  yield  K.  - 0.091,  K,  “ 0.^09, 

K.  - 1.818,  K.  - 0.182,  K,  - 2.000  and  K - 1.818. 
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In  Case  11  the  curves  tor  o and  o,  vs.  r at  t » 0 and  ■» 

r 8 

are  shown  on  Fig.  3 (identical  with  Case  l.a),  the  curves  for  c and 

0 ^ vs.  t at  r “ 1.0  are  given  as  dashed  lines  on  Figs.  4,5,  and 

lastly  the  graphs  of  v at  r “ 1.0  are  presented  separately  in 

Fig.  11.  Figures  3,4  and  5 clearly  demonstrate  that  the  asymptotic 

values  of  o and  at  the  interface  are  identical  in  Cases  l.a 

r 0 

and  II,  which  is  as  expected  since  Case  11  reduces  to  l.a  as  the  tran- 
sient effects  die  out.  We  also  note  in  Figures  4 and  5 that  the  rate 
of  approach  to  the  asymptotic  stress  values  is  generally  slower  in 
Case  II  than  in  Case  l.a,  although  in  the  special  case  m - 0 we 
again  have  t - 1 . On  comparing  the  velocity  solutions  for  these 
two  cases  (Figures  b and  11),  we  see  that,  while  vtl,t)  is  constant 
for  m = 0 and  slightly  decreasing  for  m » 2 in  Case  l.a, all  velocity 
solutions  decrease  through  a much  greater  range  in  Case  II.  This  is 
a direct  consequence  of  the  fact  that  the  two  cases  have  different 
initial  velocities  (Case  II  much  greater)  but  equal  asymptotic  veloci- 
ties, as  discussed  previously  in  section  3.2. 

We  finally  turn  to  Case  III,  in  which  the  transient  creep  is 
ignored  while  the  creep  properties  are  allowed  to  vary  with  the  radial 
position.  You  will  recall  that  we  simplified  the  analysis  in  section 
3,3  by  assuming  incompressibility  (v  “ 0.5)  and  by  employing  a multi- 
ring approximation.  In  a fuel  rod  the  temperature  and  flux  and 
accordingly  the  creep  parameters  V and  T will  decrease  in  the 
radial  direction.  For  the  purpose  of  illustration  we  shall  assume 


that  I * Td  “ T where  T is  a parabolic  function  [i.e.,  M = 2 in 
1 K 

eqs.  (23)]  with  dl'/dr  = 0 at  r = .3  , with  r(.3)/T(l)  * 2.6,  and 

fl 

with  the  mean  (1/.7)  I'(r')dr'  = 1 as  previously  required. 

J .3 

We  thereby  obtain 

r(r)  *>  1.11b  + .948r  - 1.580r2  (101) 

whereby  T(.3)  = 1.258  and  F(l)  * .484,  which  is  a change  of  moderate 
magnitude.  We  consider  that  the  average  material  behavior  is  the  same 
as  in  Case  l.a  [see  eqs.  (98)],  and  employ  the  two  ring  approximation 
for  the  radial  dependence,  with  T the  average  value  in  the  inner 
ring  (0.3  r S.  0.65)  and  r~  the  average  value  in  the  outer  ring 
(0.65  <.  r <_  1.0).  The  complete  set  of  nondimens ional  creep  constants 
for  use  in  eqs.  (93)  is  then  given  as  = 0.091,  K,  = 0.909, 
r:  = 1.193,  r7  = 0.806. 

For  Case  III  we  evaluated  or  by  means  of  the  simple  estimate 
given  by  Eq.  (4)  (which  uses  the  step  approximation  to  A(r)  through- 
out) and  v by  means  of  eq . (97);  the  curves  for  c ( r,“),  or(l,t)  and 
v(l,t)  (for  m = 0,2)  are  shown  as  dash-dot  lines  on  Figures  3,4  and 
6 respectively.  Note  that  as  expected  the  slope  of  the  radial  stress 
is  discontinuous  in  Figure  3 at  the  point  r = 0.65.  On  comparing 
Case  III  with  Case  I. a we  see  from  Figure  4 that  the  asymptotic  value 
of  f is  greater  in  Case  III,  and  thus  for  m = 0 the  contact 
pressure  in  Case  III  is  no  longer  equal  to  unity.  Finally,  on  examining 
Figure  6 we  see  that  neither  the  initial  nor  the  asymptotic  velocities 
are  equal  for  the  two  cases,  and  hence  for  m = 0 the  velocity  at  the 
interface  is  again  no  longer  constant  for  Case  111.  We  must  point  out 
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that  these  curves  are  solely  for  the  purpose  of  illustration,  and  that  it 
would  generally  be  necessary  to  use  more  than  two  rings  to  obtain  good 
estimates  for  the  contact  pressure  and  the  integration  function  g0(t)  . 
Furthermore,  as  we  have  pointed  out,  one  could  not  only  obtain  a better 
estimate  for  the  radial  stress  but  could  also  obtain  a reasonable  continu- 
ous estimate  for  the  tangential  stress  by  using  the  actual  radial  distribu- 
tion for  the  creep  properties  in  eqs.  (63). 

As  previously  stated,  our  major  goal  here  was  to  obtain  some  analyti- 
cal solutions  for  various  special  cases,  which  could  then  be  compared  with 
numerical  results  from  existing  computer  codes.  In  order  to  obtain  such 
solutions  for  this  very  complex  problem,  we  found  it  necessary  to  make  a 
number  of  simplifying  assumptions.  As  one  attempts  to  relax  these  assump- 
tions, the  complexity  of  the  problem  increases  only  slightly  in  some  cases 
while  considerably  more  in  some  other  cases.  For  example,  the  solutions 
presented  may  be  easily  extended  to  include  time  dependence  in  the  inner 
radius  and  pressure,  and  an  external  pressure  may  also  be  easily  included. 

On  the  other  hand,  the  inclusion  of  swelling  and  thermal  expansion  would 
result  in  a substantial  increase  in  complexity  since  nonlinearity  precludes 
the  use  of  simple  superposition.  Finally,  the  relaxation  of  some  other 
assumptions  would  result  in  problems  which  are  orders  of  magnitude  more 
difficult.  In  this  category  we  have  the  inclusion  of  nonlinearity  in  the 
thermal  creep  of  the  fuel  as  well  as  of  the  cladding,  since  the  compatibility 
relation  in  the  stress  function  is  then  a nonlinear  partial  differential 
equation. 
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Figure  1 - Geometry  of  the  Fuel  Rod 
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Figure  3-  Stresses  (crr  and  crd)  vs  r- Cases  l a,  H,and  HI 
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Figure  5 -Tangential  Stresses  vs. I at  r =1.0 -Cases  la  and II 


Figure  6 -Radial  Velocity  v vs.  tat  r=1.0~  Cases  la  and  HI 
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Figure  7-  Stresses  (crr  and  o^)  vs.  r-Case  l b 
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Figure  10-  Radial  Velocity  v vs. tat  r = 1.0 -Care  l b 
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